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ABSTRACT 

Context. Dust emission and scattering provide tools for the study of interstellar clouds. Detailed knowledge of grain properties is 
essential for correct analysis of these observations. However, dust properties are suspected to change from diffuse to dense regions. 
The changes could affect our interpretation of observations from large scales (e.g., clump mass spectra) down to the structure of 
individual cores. 

Aims. We study a northern part of the Corona Australis molecular cloud that consists of a filament and a dense sub-millimetre core 
inside the filament. Our aim is to measure dust temperature and sub-mm emissivity within the region. We also look for confirmation 
that near-infrared (NIR) surface brightness can be used to study the structure of even very dense clouds. 

Methods. We extend our previous NIR mapping south of the filament. The dust colour temperatures are estimated using Spitzer 
160pm and APEX/Laboca 870/jm maps. The column densities derived based on the reddening of background stars, NIR surface 
brightness, and thermal sub-mm dust emission are compared. A three dimensional toy model of the filament is used to study the effect 
of anisotropic illumination on near-infrared surface brightness and the reliability of dust temperature determination. 
Results. Relative to visual extinction, the estimated emissivity at 870 jim is Kgvo = (1.3 + 0.4) x 10~ 5 mag -1 . This is similar to the 
values found in diffuse medium. A significant increase in the sub-millimetre emissivity seems to be excluded. In spite of saturation, 
NIR surface brightness was able to accurately pinpoint, and better than measurements of the colour excesses of background stars, the 
exact location of the column density maximum. Both near- and far-infrared data show that the intensity of the radiation field is higher 
south of the filament. 

Key words. ISM: Structure - ISM: Clouds - Submillimeter - Infrared: ISM - dust, extinction - Scattering - Techniques: photometric 



1. Introduction 

The column densities of dense molecular clouds can be esti- 
mated in several ways. Because the dust can usually assumed 
to be well mixed with the gas, dust can be used as a tracer of 
cloud structure. The available tools are based on dust extinction, 
thermal dust emission, and light scattering. 

Extinction can be traced with star counts. However, in the 
NIR regime and for dense clouds it is more common to make 
use of the reddening that takes place when the light from back- 
ground stars travels through a dust cloud. The method is useful 
particularly in NIR because, at those wavelengths, the variation 
in the intrinsic colours of background stars is small, the NIR 
extinction curve is relatively insensitive to changes in the dust 
populations, and the lower optical depths allow studies of clouds 
with visual extinctions up to ~ 20™. The reliability can be fur- 
ther increased by using several NIR bands, the combination of 
J-, H-, and K-bands being the most common one (e.g., Lombardi 
& Alves|2002 Lombardi l2009l . 

A complementary tool is provided by thermal dust emission 
at far-IR and millimetre wavelengths. The emission is optically 
thin and, therefore, capable of probing higher densities than what 
are accessible by observations of background stars. The conver- 
sion between intensity and column density requires knowledge 
of dust temperature and dust emissivity. The colour temperature 
can be estimated from observations of two or more wavelengths. 
However, there are always temperature variations along the line- 
of-sight and, because warm dust produces much stronger signal 



than cold dust, the derived colour temperature can be signifi- 
cantly higher than the actual temperature of the bulk of dust. 
The problem is acerbated in pre-stellar cores where FIR obser- 
vations can completely ignore the cold core where the tempera- 
ture can be below 10 K or, at the very centre, even close to 6 K 
(Crapsi et al. 120071 Harju et al. 120081 . In geometrically simple 
cases, modelling can be used to estimate the temperature vari- 
ations and to correct the column density estimates accordingly 
(e.g., Stamatellos et al. 12007] ! 

Dust emissivity and the spectral index appear to vary from 
region to region (Cambresy et al. 2001 del Burgo et al. 2003 
Dupac et al. 2003; Krame r et al. 120031 St epnik et al. 120031 
Lehtinen et al. 120041 120071 Ridderstad et al. 120061 . The varia- 
tions are attributed to either grain growth (Ossenkopf & Henning 
1994 Krugel & Siebenmorgen |1994l l or physical changes in the 
grain material itself (e.g., Mennella et al. 1998 Boudet et al. 
2005 ). Like temperature variations, changes in grain emissivity 
should be correlated with density and, therefore, might cause 
systematic errors in, for example, core mass spectra and the es- 
timated structure of pre-stellar cores. 

With deep NIR observations one can detect the surface 
brightness that is caused when dust grains scatter photons orig- 
inating in the general interstellar radiation field. Both numerical 
modelling and first observational studies suggest that also the 
phenomenon can be used for quantitative study of cloud struc- 
ture (P adoan et al. |2006l Juvela et al. |2006ll2008l Nakajima et al. 
2008 ). The method works best when visual extinction is above a 
few magnitudes but below ~ 20 m . The lower limit is set by what 
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is observationally feasible and the upper limit by the saturation 
of the NIR signal at high optical depths. In the absence of strong 
local radiation sources, the reliability of the method appears to 
be comparable to the other methods listed above. However, be- 
cause the surface brightness is observed in NIR, the obtained 
resolution is higher, potentially even of the order of one arc sec- 
ond. 

We continue the study of the northern filament of the Corona 
Australis cloud. In previous paper Juvela et al. (2008 ), we car- 
ried out comparison between two column density estimators, one 
based on observations of background stars and one on the scat- 
tered light. In this paper, we present new observations that ex- 
tend the previous NIR mapping south of the Corona Australis 
filament. We have also mapped the region at 870 /mi using the 
Laboca instrument on the APEX telescope. The morphology of 
the sub-millimeter map is compared with FIR and NIR obser- 
vations. With the aid of Spitzer 160 /mi data, we estimate the 
dust emissivity at 870 /mi. A three-dimensional radiative trans- 
fer model is used to study the effects of asymmetric illumination 
and the differences that may exist between the colour tempera- 
ture estimates and the true dust temperature. 



2. Observations 

2.1. Near-infrared observations 

A composite image of the near-infrared J-, H-, and Ks-band data 
is shown in Fig. [T] The northern part of the field, some 4' X 8' 
in size, was observed in August 2006 using the SOFI instrument 
on the NTT telescope. Those data cover the densest part of the 
filament and, in the north, extend to an area where the estimated 
visual extinction is ~ l m . These observations were described in 
Juvela et al.|2008] 

Further SOFI observations were made in June 2007. The new 
data consist of two, partially overlapping ~ 4' x 4' fields. These 
extend the original mapping south of the centre of the filament. 
As previously, the observations were carried out as ON-OFF 
measurements so that the faint surface brightness could be re- 
covered. Three OFF fields were selected based on IRAS images 
and 2MASS extinction maps from regions of low cirrus contam- 
ination. In order to minimize the effect of sky variations, obser- 
vations of ON- and OFF-fields were interleaved in a sequence of 
OFF - ON-2 - ON-1 . The OFF-field was varied between obser- 
vation blocks to average out any faint background gradients that 
might exist in the OFF-fields. Each observed frame was an av- 
erage of six 5 second exposures, corresponding to an integration 
time of 30 seconds. After each sequence, the telescope returned 
to the OFF-position, and, at the end of the observation block, 
took six 5 second exposures of the OFF-field. The total integra- 
tion times for each ON-field was 26 minutes in J, 37.5 minutes in 
H, and 97.5 minutes in Ks band (see Table[TJ. The observations 
were calibrated to 2MASS scale by using photometry of ten se- 
lected stars in each individual frame. Their average was used to 
correct the observed fluxes to above the atmosphere. The two 
ON-fields were mosaiced together using the overlapping area. 
To have the same signal in the overlapping area, this required 
an additional small correction to ON-1 frame. The overlapping 
area between the 2006 and 2007 observations, slightly more than 
1' x 1' in size, was used for the final adjustment of the new sur- 
face brightness measurements. 



Table 1. Positions and total integration times of the observed ON 
and OFF fields. The last column gives list for the OFF fields the 
extinction estimates of Schlegel, Finkbeiner & Davis ( 1998). 



Field 


centre position 


tQ) 


f(H) 


t(Ks) 






(J2000) 




(min) 


(min) 


(min) 


(mag; 


ON-1 


19 h m 42\ 


-37°0'54" 


26 


37.5 


97.5 




ON-2 


19 h m 47\ 


-36°57'54" 


25.5 


37.5 


97.5 




OFF-1 


18 h 58 m 14 s 


, -36°50'11" 


11 


16 


35 


0.45 


OFF-2 


18 h 58 m 15 s 


, -36°34'43" 


12 


16 


35 


0.35 


OFF-3 


18 h 56 ra 59 s 


, -36°53'5" 


5.5 


8 


35 


0.33 




19 01 19 00 3(3-0 

RA (2000,0) 



Fig. 1. A composite of the near-infrared observations of the 
Corona Australis filament. The J-, H-, and Ks-bands are coded, 
respectively, in blue, green, and red. The filament runs horizon- 
tally across the image and, because of the stronger saturation of 
the J- and H-band surface brightness, shows up as a red region. 



2.2. Sub-millimeter and far-infrared data 

The Corona Australis filament was mapped with APEX/Laboca 
at 870 /im. The first observations were carried out in November 
2007. However, the weather conditions were not optimal and, 
additionally, some detector scans were affected by the bright R 
Corona Australis region that is located some 15' west of our 
NIR field. In the following, we rely mostly on data gathered in 
October 2008. Observations consisted of scans ~30' in length. 
These were made mainly in the North-South direction, i.e., per- 
pendicular to the filament. To improve sky noise suppression, 
the position angle of the scans was varied by up to ~30 de- 
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Fig. 2. The 870 /mi Laboca map. The white line denotes the out- 
line of the NIR map (see Fig.[TJ. The bright point source T Tauri 
star S Cr A (Jov fI945l Carmona et al. 120071 



grees. The precipitable water vapour content of the atmosphere 
was ~ 1.2 mm. The observations were reduced using the Boa 
program, version 1 . lQ. The atmospheric opacities derived from 
skydips using Boa are known to be underestimated. Following 
instructions on the official Laboca web-site, the opacities were 
scaled up by a factor 1.3. 

During observations, Uranus, Neptune, and, as a secondary 
calibrator, IRAS 16293, were observed at regular intervals. The 
final maps contain one bright millimeter point source, S Cr A, 
which gives possibility for another check of calibration. The 
source is a T Tauri star that corresponds to the submillime- 
ter source MMS-7 observed by Nutter et al. (2005) with the 
SCUBA instrument. At 850 /mi, their estimate of the source flux 
was 0.70+0.02 Jy/beam. Taking into account the difference in 
the beam size and in the wavelength, our estimate is higher by 
slightly less than 8%. The 870//m map is shown in Fig. [2] 

The area has been mapped in Spitzer guaranteed time pro- 
grams, at 3.6, 4.5, 5.8, and 8.0 /mi with the IRAC instrument and 
at 24, 70, and 160//m with the MIPS instrument. In this paper, 
we use only the longest wavelength, 160 /mi, to derive tempera- 
tures of the large dust grains that are responsible for the sub-mm 
signal. The shorter Spitzer wavelengths probe small grain and 
PAH populations. These should be included in a comprehensive 
model of the cloud but are not used here, because they are largely 
independent of the large grain population and, compared with 
the Laboca data, trace outer cloud layers. 



3. Results 

3.1. NIR observations 



In paper I (Juvela et al. 2008 ), good agreement was found be- 
tween Ay estimates obtained from the reddening of background 
stars and from the surface brightness of scattered light. With the 
new observations, the comparison can be extended to the south- 
ern side of the filament. For the background stars, Ay estimates 
were derived using the NICER method (Lombdardi & Alves. 
I20011 I. The surface brightness was converted to Ay using the 
same method and exactly the same parameters as in paper I. 
The extinction maps are shown in Fig. [3] while Figure [4] shows 
North-South profiles of the filament in Ay and scattered surface 
brightness. 



see http://www.astro.uni-bonn.de/boawiki 



As concluded in paper I, on the northern side (offsets 0- 
5') the methods give very similar extinction values. Within 
the filament the surface brightness becomes gradually saturated 
(see Fig. [TJ and the uncertainty of the Ay estimates increases. 
Therefore, for the surface brightness data, predictions are shown 
only up to Ay = 15™. Within the filament, the density of back- 
ground stars decreases. This leads to bias also in the NICER esti- 
mates, as discussed in Paper I. Furthermore, in the region where 
the Ks-band surface brightness is completely saturated also no 
background stars are seen. Therefore, the true peak extinction 
can be higher than shown by the NICER values. 

On the southern side, the surface brightness predicts ~20- 
30% higher extinctions than the colour excesses of the back- 
ground stars. Compared to the NICER extinction profile, espe- 
cially in the J- and H-bands, the surface brightness is higher 
on the southern side of the filament (see Fig. [4] lower frame). 
This suggests that the intensity of the radiation field would be 
higher on the southern side of the filament. The conversion be- 
tween scattered intensity and Ay depends on assumptions of dust 
properties and the illuminating radiation field. In Paper I, the pa- 
rameters were fixed using data north of the filament. Therefore, 
assuming identical dust properties, the difference in the Ay esti- 
mates would indicate a similar percentual difference in the inten- 
sity of the NIR radiation field. In the K-band the asymmetry is 
smaller than in the other two bands. This suggests a change in the 
spectrum of the illuminating radiation although the NIR colour 
ratio might also be affected by observational uncertainties. 

Both the new NIR observations and the existing far-infrared 
data (see Sect. \3.2\ indicate that the intensity of the radiation 
field increases towards south. In paper I an opposite conclu- 
sion was reached on the basis of NIR colours close to the dip in 
the Ks-band intensity. The intensity ratio J/Ks was found to be 
higher in the northern part of the dip which, because the cloud is 
more transparent for Ks-band radiation, would seem to be con- 
sistent with most of the illumination coming from that direc- 
tion. However, the situation may be more complicated and the 
outcome depends on whether the observed local minimum of 
K-band surface brightness exactly coincides with the extinction 
maximum or is shifted away from the direction of stronger illu- 
mination. We return to this question in Sect.|4]where we examine 
the situation with the help of a radiative transfer model. 

3.2. Comparison of NIR, far-infrared, and sub-mm emission 

Figure [5] compares the morphology of Cr A filament as seen at 
different wavelengths. The colour image shows the Ks-band sur- 
face brightness that should be well correlated with the true col- 
umn density up to Ay ~ 10 m but which, within the filament, 
becomes saturated and exhibits a dip at the centre of the fila- 
ment. In addition to the NIR data, the figure shows Spitzer 160^ 
observations and our Laboca observations at 870 /mi. 

The sub-mm dust emission shows very good correlation with 
the morphology of NIR emission. At the location of 870 /mi 
maximum, the NIR signal is saturated but, based on the ratio 
of J, H, and Ks intensities, this was already concluded to be the 
position of the highest column density (Juvela et al. 2008]). In 
other words, the dip in the Ks-band image pinpointed the loca- 
tion where the sub-mm emission was going to peak. 

In Sect. |3.1| we concluded that the radiation field is stronger 
on the southern side of the filament. In Fig.|5]the 870pm emis- 
sion is displaced by some 10" south-east of the Ks-band dip. 
This is consistent with the idea of asymmetric illumination: if 
the dip in NIR surface brightness is shifted away from the di- 
rection of stronger radiation, the sub-mm peak should similarly 
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HA f 20O0.0) RA 

Fig. 3. Extinction maps obtained using the background stars (the NICER method; frame a) and observations of scattered light 
(frame b). The black contours are drawn at 2, 4, 6, 8, 10, and 15 and the white contours at 20 and 30 magnitudes of visual extinction. 
The resolution is 20" for the NICER map and 10" in for the extinction maps derived from surface brightness measurements. 



be shifted slightly towards the stronger radiation. In Fig. [5] the 
asymmetry points towards the nearby T Tauri star S CrA as 
the cause of these shifts. However, at 160 /im the emission is 
very extended on the souther side of the filament. This reflects 
also the mass distribution that is asymmetric with respect to the 
densest part of the filament. In the 160yimi emission, the peak is 
directly south of the sub-mm peak, not towards the star S Cr A. 
At shorter wavelengths, the 100 //m surface brightness rises from 
~20 MJy sr~' in the northern part of our NIR field to 41 MJy sr -1 
in the filament. However, south of the dense filament, the surface 
brightness stays above 35 MJy sr -1 . In the Spitzer 70 /mi maps, 
the surface brightness rises as one moves southwards across the 
filament, without clear correlation with the position of S Cr A. 

3.3. Dust colour temperature in the Corona Australis filament 

The colour temperature of large grains was estimated using the 
160 yum Spitzer MIPS map and our APEX/Laboca 870 /mi ob- 
servations. The values are based on the assumption of a single 
gray body spectrum, B y (T) x . Unfortunately, the value of the 
spectral index is uncertain (e.g., Dupac et al. 2003) and, because 
dust properties may exhibit variations that are correlated with 
density and temperature, /3 may have significant variations even 
within a single cloud. There is also indication that the spectral 



index is wavelength dependent and that, compared to sub-mm 
and millimetre wavelengths, the variation of /3 is stronger in the 
FIR regime. With the available data, we cannot study these vari- 
ations in detail. The colour temperature map obtained for a fixed 
value of ft = 2.0 is shown in Fig. [6] The figure also indicates 
the area that was used for the determination of the intensity zero 
level at the two wavelengths. 

In the reduced 870 /mi map, only the densest filament is vis- 
ible. On both sides, the intensity first becomes slightly negative 
before again raising towards north and south. This suggests that, 
in the data reduction, the baseline subtraction was not perfect. 
However, the overshoot to negative values is only ~5% of the 
peak intensity and, for the colour temperature calculation, the 
zero level is determined from an area where the surface bright- 
ness reaches minimum (see Fig. [6]). Therefore, the uncertainty 
in the baseline subtraction translates to an uncertainty in colour 
temperature that is only ~0. 1 K for the sub-mm peak but can be 
close to one degree in areas of lower intensity where the colour 
temperature is close to 15 K. A related source of uncertainty re- 
sults from the fact that, in the 870 /mi map, some of the large 
scale structure may have been filtered out. 

In the filament, the colour temperature is reduced to 
Tc =10.7 K. In this value, the uncertain value of /3 is a larger 
source of uncertainty than the observational errors. Varying f3 in 
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Fig. 4. Upper frame: Comparison of Ay profiles of the cloud fil- 
ament derived using the NICER method and using the scattered 
surface brightness. For surface brightness, values are truncated 
at Ay = 15 m because of the saturation of the scattered signal. 
Lower frame: The profile of the filament in three NIR bands. The 
cut is from -36°49'12" in the north (offset zero) to -37°2'21" in 
the south. Values are averaged over the interval from 19 h m 43.8" 
to 19 h m 50.4". 




B* dJOOO.O) 

Fig. 5. The Ks-band surface brightness (colour image) with 
contours of 870/itn (white lines) and 160/vm (black lines) sur- 
face brightness. For 870jum, the contours are from 0.02 to 
0.12mJy/beam, at intervals of 0.02mJy/beam. The 160 //m con- 
tours are drawn at 140, 160, 180, and 200 MJy sr -1 . 
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Fig. 6. The 870 pm Laboca map. The white line denotes the out- 
line of the NIR maps (see Fig . [TJ . The bright point source T Tauri 
star SCr A (Jov [19451 



the plausible range from 2.5 to 1.5 would change the minimum 
colour temperature between 9.5 and 12.4 K. A change AB-0.5 
corresponds to a change by a factor of 2.7 in the relative dust 
opacity at 160 ^m and 870 pm. Because of temperature varia- 
tions along the line-of-sight, the colour temperature is likely to 
overestimate the mass-weighted true dust temperature. This bias 
is examined further in Sect. [4] in connection with a radiative 
transfer model. 



3.4. Sub-mm dust emissivity 

Using the derived colour temperatures Tc and the surface bright- 
ness I v , the dust emissivity at 870 pm, relative to visual extinc- 
tion, can be calculated differentially, 



AI v /B v (Tc) . _u 

Kfno = — ; mag 

AA V 



(1) 



This assumes that dust emissivity remains constant over the 
whole area. Figure [TJa shows the correlation between I v /B v (Tc 
and visual extinction. To lower the noise of Ay, we use weighted 
average of the maps shown Fig. [3] At low Ay, equal weight is 
given to the two maps but, by Ay ~ 23™, the values are based 
on NICER method alone. The slope of the linear fit gives a value 
#870 - 0.6 x 10~ 5 mag _1 . However, the slope steepens towards 
the filament and, for Ay > 13 m , one obtains emissivity that is 
twice as high, <f870 = 1.3 x 10~ 5 mag -1 . 

Emissivity can also be estimated directly on the basis of 
absolute values. As above, the values are low outside the sub- 
millimeter peak but increase to ~ l.lxlO" 5 mag -1 at the location 
of the surface brightness maximum. In this case, the zero level 
of surface brightness values was estimated north of the filament, 
around declination -36°53'. 

There is a temptation to interpret the spatial variations as an 
indication of increase of dust sub-mm emissivity towards the 
densest parts of filament. However, the 870 pm surface bright- 
ness is detected only close to the actual sub-mm peak and the 
emissivities, obtained either from correlation or absolute surface 
brightness, are uncertain outside this region. Furthermore, in the 
870 pm map, some of the large scale structure may have been fil- 
tered out. This could account for the lower value obtained given 
by the absolute surface brightness values and could contribute 
to the decrease seen outside the surface brightness maximum. 
Therefore, considering both the S/N ratio and the possible filter- 
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Fig. 7. Ratio of 870 yum surface brightness and the Planck func- 
tion, B(Tc), as function of visual extinction. Least square fits 
are shown for all data with 7 nll (870) > 0.01 Jy/beam and for the 
central part of the filament where Ay > 13 m . 



ing out of large scales, the values obtained at the very centre of 
the filament should be the most reliable ones. 

The emissivities were calculated using the 160^/m/870/im 
colour temperature which may overestimate the true dust tem- 
perature if there are strong temperature variations along the line- 
of-sight. In Sect. |4] we will examine, whether this could lead to 
underestimation of k&jq. In paper I, we concluded that, extinction 
within the filament may be underestimated because of the steep 
extinction gradients. The bias, which was estimated to be ~25% 
at Ay ~ 20™, would cause k&jq to be overestimated, probably by 
a similar factor. We quote emissivities relative to visual extinc- 
tion. In the NICER method, the conversion to Ay was done as- 
suming a standard extinction curve that corresponds to Ry=3.1. 
Although it is probable that Ry increases towards the dense fila- 
ment, the effect on K/Ay is only of the order of one percent, i.e., 
negligible compared with the other sources of uncertainty. 

Because of the various sources of uncertainty, we can only 
conclude that, at the centre of the Cr A filament, the dust emis- 
sivity is likely to be at least ~ 1 x 10~ 5 mag -1 but could easily 
be higher by ~5Q%. Relying on the values obtained at the cen- 
tre of the filament and taking into account the various biases, we 
write our estimate as /cg7o = (1.3 + 0.4)x 10~ 5 mag -1 . The quoted 
uncertainty is, of course, only approximative. 



4. Toy model of the filament 

A toy model of the Cr A filament was used to examine two ques- 
tions, the effect of asymmetric illumination on the observed NIR 
surface brightness and the bias of the dust colour temperature 
used in the derivation of sub-mm dust emissivity. The model 
consists of a cylinder, representing the large-scale filament, and 
a denser spherical region in the centre, representing the observed 
sub-mm peak. The density follows gaussian profile, both with re- 
spect to distance from the axis of the cylinder and from the cen- 
tre of the spherical component. The extinction is ~10 m across 
the cylinder and reaches ~ 70 m through the centre, i.e., through 
the combined density peak of the cylinder and the sphere. The 
model was resampled onto a cartesian grid of 65 3 cells. 

The dust properties are taken from Draine (2003 ). The scat- 
tering was calculating using the tabulated phase function avail- 



able on the web Q, assuming the dust model corresponding to 
Ry = 5.5. For radiative transfer calculations, we use our Monte 
Carlo programs (Juvela & Padoan 2003] Juvela [2005b . 

4.1. Correlation of NIR surface brightness 

In Paper I the spatial variations in the NIR intensity ratios were 
examined and several regions were identified that exhibited a 
particular combination of surface brightnesses. Using the toy 
model, we re-examine the interpretation of the regions c — d and 
f — g that were shown in Fig. 4 of that paper. 

Figure [8] shows the correlation between J- and K-band sur- 
face brightness in the toy model. The model cloud is illuminated 
by an isotropic radiation field that corresponds to the NIR ISRF 
given in Lehtinen et al. (|1996l l. Additionally, in the south (as 
marked in the figure) there is an additional radiation source with 
similar spectrum that raises the total intensity by 30%. This in- 
troduces asymmetry that is clearly visible as northern and south- 
ern branches in the correlation plotted in Fig. [8};. 

Figure [8}; can now be compared qualitatively with Fig. 4 of 
paper I. The regions c and d are clearly consistent with asym- 
metric illumination with stronger radiation in the south. More 
specifically, the asymmetry suggests that the additional radiation 
is coming from south-east. 

For regions / and g the correspondence with the toy model is 
less clear. In the model, the northern and southern branches (as 
marked in Fig. [8};) cross each other before the location of peak 
extinction is reached. Therefore, the intensity ratio J/K should be 
higher on that side of the Ay peak where the radiation is stronger. 
If regions / and g are interpreted as the tip of the two branches, 
they would perfectly match a case where radiation comes mainly 
from the north rather than from the south. Because other evi- 
dence clearly shows that this is not the case, the final interpre- 
tation of the NIR colours at the centre of the filament remains 
open. Smaller scale density structures (e.g., structures extending 
in front of the main cloud) could provide one possible explana- 
tion. 

In the model the southern radiation source extended 10° from 
the southern direction. Figure [8]remains qualitative similar even 
if the opening angle is as large as 45°. Similarly, a clear loop that 
is formed in Fig. [8};) when the southern and northern branches 
cross, persists even when the radiation source is moved +45° 
from the south towards the front or back of the cloud. However, 
the the separation of the two branches is much smaller if the 
source is moved towards the front of the cloud. If the source is 
moved towards the back of the cloud, the peak surface brightness 
increases but the intensity at the centre of the cloud decreases, 
both by up to ~20%. These effects suggest that, with more de- 
tailed modelling of observations, further constraints may still be 
derived for the radiation field. 

4.2. Colour temperatures from the toy model 

The determination of sub-mm emissivity can be affected by er- 
rors in the estimated dust temperature. In particular, when there 
are temperature variations along the line-of-sight, the observed 
colour temperature is known to be biased towards higher tem- 
peratures (e.g. Lehtinen et al. I2007I I. We use the toy model to 
estimate the magnitude of this bias. We use exactly the same 
model as above, i.e., including the additional radiation source in 
the south that increases the total illumination by 30%. 



2 http://www.astro.princeton.edu/~draine/dust/ 
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Fig. 8. NIR surface brightness in our toy model, a. Image of the 
K-band surface brightness. The cylindrical filaments runs verti- 
cally across the image and an additional density peak is located 
at the centre of the image. The illumination is stronger on the 
right hand side of the filament, b. Surface brightness and visual 
extinction profiles of the filament averaged over the area lim- 
ited by the dashed lines drawn in frame a.. The lower dotted 
line shows Ay cross section of the cylinder alone, c. Correlation 
of J- and /iT-band surface brightness for pixels between dashed 
lines in frame a. Colours highlight pixels close to the Ay peak. 
Yellow triangles (tip pointing upwards) correspond to pixels on 
the side of higher illumination and cyan triangles (tip pointing 
downwards) to pixels away from the stronger illumination. In 
frames a and c, the letters N and S refer to north and south, fol- 
lowing the hypothesis that the illumination of the Cr A filament 
is stronger in the south. 

Figure [9}3 shows a map of colour temperatures that is calcu- 
lated from modelled surface brightness at 160 /mi and 870 fim. 
A value of B = 2.0 is used for the emissivity index. This is very 
close to the spectral index of the employed dust model, which is 
1.96 between the two frequencies. Frame b shows the difference 
between the observed colour temperature and the mass-weighted 
average temperature along each line-of-sight. As expected, the 
colour temperature overestimates the true mean temperature of 
dust grains. The error increases with opacity and is just over 2 K 
at the centre of the cloud. 

Because the dust emission is optically thin, the ratio of the 
true column density and the column density estimated on the 
basis of the colour temperature can be calculated from 



Jo 



(s)B v (T(s))ds/[B v (T c ) n(s)ds] 



f n( 
Jo 



(2) 



Here T is dust temperature, Tq colour temperature, and s dis- 
tance along the line-of-sight. Figure|9]; shows this ratio as func- 
tion of observed colour temperature. At the cloud centre, the true 
column density is more than ~5Q% higher than what is estimated 
on the basis of colour temperature. This suggests that the value 
of K87o derived in Sect |3.4| could be underestimated although not 
necessarily by the same amount. 

To estimate the effect on k^q, we scaled the 870//m surface 
brightness values upwards according to the relation of Fig. [9}:. 
If the CrA filament behaves as our model, the scaled surface 
brightness should correspond to a situation where all dust along 
a line-of-sight has temperature equal to the colour temperature. 
Repeating the analysis of Sect. 13.41 values of k^q increase less 
than 10% in the fit to all data and ~20% in the fit to data above 
Ay = 13 m . The increase in the peak value of Fig. is similarly 




Fig. 9. Comparison of observed colour temperatures and the true 
dust temperature in the toy model, a. The dust colour tempera- 
ture based on 160 //m and 870 fim surface brightness, b. The dif- 
ference between the colour temperature and the mass weighted 
average temperature along each line-of-sight. c. The ratio of 
true column density and the column density estimated based on 
the observed colour temperature. The plot shows pixels that, in 
frame b, fall between the dashed lines. Pixels south of the sym- 
metry axis are plotted as red plus signs. 



~20%. While noticeable, the bias caused by the use of colour 
temperatures may be compensated by the opposite bias in Ay 
values. Therefore, the value derived in Sect. 13.41 ^70 = (1.3 ± 
0.4) x 10~ 5 mag -1 , remains our best estimates. 



5. Discussion 

The IRAS and Spitzer far-infrared observations demonstrate an 
asymmetry where the intensity of the radiation field and/or the 
extent of diffuse material is larger south of the Corona Australis 
filament. Our extinction measurements show that the extinction 
is reduced below Ay ~ 3 m within 5' of the filament, both in 
the North and in the South. The column densities are slightly 
larger in the south but NIR and far-infrared colours indicate 
that the asymmetry is mainly caused by the illumination of the 
cloud. Our toy model (Fig. [HJ showed that, at large scales, the 
NIR colour variations are consistent with this interpretation. In 
the model, the radiation field was increased by 30% by adding 
one radiation source in the south but this resulted in larger NIR 
colour variations than what was observed. Either the asymmetry 
is, at NIR wavelengths, less strong or the radiation comes from 
a much large solid angle or from a direction that deviates signif- 
icantly from the exact southern direction. 

Stars S CrA and HD 176269 lie a few arc minutes south- 
east of the filament. HD 176269 is a pre-main sequence B9V 
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star which, according to the HIPPARCOS parallax, is at the dis- 
tance of Corona Australis, ~ 140 pc. In Digital Sky Survey Blue 
map, there is a reflection nebula around the star. The 2MASS 
magnitudes of the star are 6.7, 6.8 and 6.7™ for the J, H and Ks 
bands, respectively. Assuming the intervening extinction is simi- 
lar to the extinction between Cr A and us and taking into account 
the projected distance between the star and the outskirts of the 
filament, ~ 3', the star could produce an additional intensity that 
is ~ 3.4(2.1, 1.4) x 1(T 18 ergs" 1 cm" 2 Hz" 1 in J-band (H-band, 
Ks-band). This would add a significant, bluer component to the 
ISRF. 

The 70 /im surface brightness shows a clear surface bright- 
ness gradient across the filament and brightening towards the 
south. This supports the assumption that the additional illumi- 
nation comes directly from the south, rather than in more to- 
wards the front or back of the filament. In the south, the ratio of 
the J- and H-bands to the Ks-band decrease from east to west, 
this also being consistent with bluer illumination from the star 
HD 176269. S Cr A is located much closer to the dense part of 
the filament, but its spectrum is redder and would not explain 
the different colour of scattered light further south of the fila- 
ment. The IRIS 100 jum map implies that the influence of S Cr A 
to the heating of dust is restricted to its immediate surroundings 
but it could have a contribution to the NIR flux impinging on the 
filament. 

In paper I we correlated NIR surface brightness with the 
extinction derived on the basis of the colour excesses of back- 
ground stars. We concluded that the diffuse signal consists 
mainly of scattered light and it can be used for quantitative study 
of cloud structure. Unless the method is calibrated using colour 
excess data (see Juvela et al. 2008; Nakajima et al. 120081 1, the 
conversion to column density requires assumptions of the radi- 
ation field and dust properties. The method should remain reli- 
able as long as the radiation field is roughly constant within the 
examined region. In the Cr A field, this condition is clearly not 
fulfilled because of the high opacity of the central filament and 
the asymmetry of the illuminating radiation. In the south, assum- 
ing similar conditions as on the northern side of the filament, the 
column density derived from NIR surface brightness is overes- 
timated by ~20% (see Fig. |4^). However, the results from the 
surface brightness data can always be corrected with the aid of 
colour excess measurements. This can be done at any scale down 
to the resolution allowed by the spatial density of background 
stars. Compared to the use of background stars only, this should 
result in a more accurate column density measurements and, as 
far as allowed by the observational noise of surface brightness 
data, extension to smaller spatial scales. 

Our estimate for the 870 /im emissivity was a^o = (1.3 ± 
0.4) x 10 -5 per magnitude of visual extinction. In the dust mod- 
els employed in our toy model, the 870 yum dust extinction cross 
sections are very similar, 1.3 x lO^mag-' for Ry = 3.1 and 
1.2 x 10~ 5 mag-' for Ry = 5.5. In this respect, the models are 
consistent with the actual filament. 

Sub-mm data obtained from PRONAOS experiment favours 
dust spectral index values close to /3 — 2 (e.g., Bernard et 
al. 119991 Stepnik et al. 2003). The spectral index is, however, 
temperature dependent and rises to ~2 only when temperatures 
is below ~20K (Dupac et al. I2003I I. There is some scatter in 
the values and, even at low temperatures, the spectral index may 
vary within a wide range from -1.5 to ~2.5. The determination 
of the spectral index is, of course, closely related to the dust 
emissivity. As mentioned in Sect. 13.31 varying the value of ft in 
the range from 1.5 to 2.5 will cause the minimum temperature 
to vary between 9.5 and 12.4 K. In the value of k, this corre- 



sponds to an uncertainty of ~30%. If sub-millimetre emissivity 
had shown significant increase, the true spectral index [5 would 
probably have been smaller than assumed. By overestimating /3, 
the dust temperature would have been underestimated and the 
obtained value of would have been even higher than the true 
value. 

The obtained value of ^70 is similar to the values that 
have been obtained in low density regions. Using COBE data, 
Boulanger et al. ([1996) derived effective dust cross section for 
dust in high latitude diffuse medium, 

o- H = 10- 25 (/l/250^m)- 2 cm- 2 . (3) 

Using a conversion factor N(H2)/Ay = 0.94 x 10 21 cirr 2 mag -1 
(Bohlin et al. 1978), equally valid for diffuse regions, this corre- 
sponds to o-(870;um) = 1.55 x 10 -5 mag-'. 

In the literature, higher emissivities are often reported for 
dense and cold clouds. Stepnik et al. (2003) studied a filament 
in Taurus. According to star counts, the extinction was at least 
A v ~ 4 m at the resolution of 3.5 '. However, using 2MASS 
H-K colours, Padoan et al. (2002) estimated the peak extinc- 
tion to be ~12 m . Based on PRONAOS observations, Stepnik et 
al. concluded that the sub-mm emissivity of big grains had in- 
creased by a factor of ~3.4 within the filament. At shorter wave- 
lengths, Lehtinen et al. (2007) reported a similar trend in the 
cloud L1642. The ratio of 200 /im optical depth and visual ex- 
tinction increased by a factor of four when moving from re- 
gions with 7dust ~ 20 K to dense regions with 7d us t ~ 12 K. 
Kramer et al. (2003) estimated variations of 850yum emissivity 
in dark cloud IC 5 146. The emissivity was found to change from 
k&5o/kv ~ 1.3 x 10- 5 in warm regions to ~ 5 x lO -5 in colder re- 
gions with T ~ 12 K. On the other hand, in their study of another 
Taurus filament, Nutter et al. (2008) found no such strong vari- 
ations. The observations, including Spitzer 70//m and 160//m 
and SCUBA 450 fim and 850 fim data, could all be modelled ac- 
curately without resorting to any spatial variations of dust prop- 
erties. The filament was similar to our Cr A field and had higher 
extinction than the region studied by Stepnik et al. 120031 

In our 870 fim observations of the Cr A filament, part of large 
scale structures may have been filtered out and this could lead 
to underestimation of dust emissivity. However, this should not 
have a large effect at the short spatial scales used in Sect. |3.4| Our 
observations suggest that the emissivity in the Corona Australis 
filament is similar to the values of diffuse medium. In spite of 
several sources of uncertainty, it is very unlikely that the emis- 
sivity could be even twice as high. This is not necessarily in 
strong contradiction with the results of Stepnik et al. (2003)) be- 
cause our K870 is an effective value for the whole line-of-sight. In 
the case of Stepnik et al., the emissivity of large grains was in- 
creased by a factor 3.4 only at the centre of their radiative trans- 
fer model. Although the observed value of k should be domi- 
nated by this densest region, the effective value for the whole 
line-of-sight should still be somewhat lower. However, the stud- 
ied CrA filament appears to be clearly different from the dense 
cores of IC 5146. In the study of Kramer et al. (2003), ksso/kv 
was above 2 x lO -5 everywhere the dust colour temperature 
(based on 450/mi and 850yimi observations) was below 14 K. 
These results indicate that the changes in sub-mm emissivity are 
not a universal function of density or temperature. 

In a future paper, we will construct a three-dimensional ra- 
diative transfer model for the Cr A filament. Varying the density 
distribution, radiation field, and dust properties, we will seek for 
a solution that would be consistent with observations from NIR 
to sub-mm. 
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6. Conclusions 

We have complemented our previous NIR observations of the 
northern filament of the Corona Australis molecular cloud. The 
dust extinction and NIR scattering are now mapped both north 
and south of the filament. The observations, together with exist- 
ing far-infrared data, show that the illuminating radiation field 
becomes stronger in the south. 

The area was mapped at 870 //m using the APEX Laboca 
instrument. With the help of Spitzer 160//m maps, we esti- 
mate dust colour temperatures which, within the central fila- 
ment, decreases down to Tc ~ 10.5 K. On the basis of the colour 
temperatures, the dust emissivity at 870 //m is estimated to be 
^870 = (1.3 ± 0.4) x 10~ 5 mag -1 as a value relative to Ay. In spite 
of the high visual extinction, the estimate is close to the values 
found in diffuse regions. 

A three-dimensional toy model of the filament was pre- 
sented. The model demonstrates that the main NIR colour 
variations described in paper I can be explained as a result 
of anisotropic illumination. The differences between observed 
colour temperatures and mass-averaged true temperatures of 
large dust grains were examined using the same model. The 
use of colour temperatures will lead to underestimation of k^q 
values but the effect is no more than ~25%. In the case of the 
Corona Australis filament, this could be compensated by an op- 
posite bias in the Ay estimates. 
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